Mean Molecular Weight 


The thermodynamic relations between P, p, and T, as well as the 
calculation of stellar opacity requires knowledge of the system’s 
mean molecular weight (defined as the mass per unit mole of 
material, or, alternatively, the mean mass of a particle in Atomic 
Mass Units). Recall that a mole of any substance contains N4 = 
6.02252 x 1078 atoms. Thus, the number density of ions is related 
to the mass density, p, by 


N 


HMa H 
where m, is the mass that is equivalent to 1 A.M.U. If the mass 
fraction of species 2 is x;, then its number density is 
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where A; is the atomic weight of the species. The number density 
of all ions in a volume of gas is then 
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To compute the contribution of (massless) electrons to the mean 
molecular weight, let Z; be the atomic number of species 7, and f; 
be the species’ ionization fraction, i.e., the fraction of electrons 
of i that are free. The number density of electrons is therefore 
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Note that in the case of total ionization (f; = 1), this equation 
simplifies greatly. Since Z;/A; = 1 for hydrogen, and ~ 1/2 for 
everything else, 
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From the definitions above, the total number density of particles 
is N 
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where the mean molecular weight is defined as 
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The Ionization Fraction 


The calculation of mean molecular weight requires knowledge 
of the chemical composition of the material and the ionization 
fraction. To calculate ionization fraction, one needs the Saha 
equation. 


In general, the Saha equation can be used to compute ionization 
fractions over most of the star. It does, however, require that 
the gas be in thermodynamic equilibrium. ‘This is true through- 
out the star, as at high densities, collisions will control the level 
populations. This approximation only breaks down in the solar 
corona, where the densities become very low. 





The Saha equation also breaks down in the centers of stars, where 
high densities cause the ionization energies of atoms to be re- 
duced. (Obviously, if the mean distance between atoms is d, then 
there can be no bound states with radii greater than ~ d/2.) In 
the case of the hydrogen atom, the Bohr radius of level n is 
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Thus, if the particle density is 
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then all the hydrogen is necessarily pressure ionized. In practice, 
the Saha equation begins to break down at nuclear distances of 
~ 10 ao, which corresponds to ~ 2.7 x 1073 u g-em~?. To correct 
for this effect, the Saha equation is normally used until it begins 
to show decreasing ionization fractions toward the center of the 
star. When this happens, complete ionization is assumed. 


To derive the Saha equation, begin by considering the Boltzmann 
equation, which states that the number of atoms in level 7 relative 
to level 7 is 
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where w; is the statistical weight of the level 7 (i.e., the number 
of separate, individual states that are degenerate in energy), and 
Xij is the difference in energy between the two levels. The number 
of atoms in level ¿ relative to the number in all levels is thus 
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where y; is the energy difference between the it” level and the 
ground state. The variable u is the partition function for the 
atom (or ion). Because u is a function of temperature, it is 
sometimes written u(T). 


Now let’s generalize this equation to electrons in the continuum. 
Let n; be the number of atoms in all levels (defined as n above), 
and let state z+1 be that where an excited electron is in the con- 
tinuum with momentum between p and p+ dp. The Boltzmann 
equation then gives 
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where y; is the energy needed to ionize the ground state of the 
atom, and dw;41 is the statistical weight of the ionized state. 


Now consider that dw; has two components: one from the ion 
(wi+1), and other from the free electron (dwe). The former is 
just the statistical weight of the ground state of the ion, while 
the latter can be computed using the exclusion rule. Since each 
quantum cell in phase space can have only two electrons in it 
(spin up and spin down), then the number of degenerate states 
in a volume h? is 
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Thus 
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The number of electrons in volume fdV = 1/ne, so the total 
number of electrons in all continuum states is therefore 
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or, if we let 2? = p?/2m_-kT, then 
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Finally, note that for the calculation above n;,1 represents those 
atoms of species n; that have one electron in the continuum state, 
i.e., ionized. It does not consider atoms of n;+ı that are them- 
selves excited. (In other words, nj41 in (5.2.4) only includes 
ionized atoms in their ground state.) To include all the excited 
states of n;41, we must again sum the contributions in exactly 
the same way as we did in (5.2.2). Thus, the statistical weight 
in (5.2.4) should be replaced by the partition function, and 
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This is the Saha equation, which relates the number of atoms in 
ionization state 2 + 1 to the number in ionization state i. Note 
that if need be, we can substitute the electron pressure for the 
electron density using P, = nekT, and write the Saha equation 
as 
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The sense of these equations is intuitive: the higher the tem- 
perature, the greater the ratio, but the higher the density (or 
pressure), the lower the ratio due to the greater possibility for 
recombinations). 





